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Recently, we have proposed a method for the local detection of quantum correlations on the basis 
of local measurements and state tomography at different instances in time [Phys. Rev. Lett. 107, 
180402 (2011)]. The method allows for the detection of quantum discord in bipartite systems when 
access is restricted to only one of the subsystems. Here, we elaborate the details of this method and 
provide applications to specific physical models. In particular, we discuss the performance of the 
scheme for generic complex systems by investigating thermal equilibrium states corresponding to 
randomly generated Hamiltonians. Moreover, we formulate an ergodicity-like hypothesis which links 
the time average to the analytically obtained average over the group of unitary operators equipped 
with the Haar measure. 
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I. INTRODUCTION 

The field of quantum information theory is dedicated 
to developing computational techniques with an advan- 
tage over classical methods using the laws of quantum 
mechanics [1] . A variety of tools for communication and 
computation science have been developed in the past 
years, ranging from quantum teleportation [2] and quan- 
tum dense coding [3] to efficient algorithms for quantum 
computers [4-6] . The fundamental resource for these ap- 
plications is usually summarized under the term quan- 
tum correlations, even though it has proven difficult to 
identify a common resource to all of these applications. 
More precisely, ideas like quantum teleportation and the 
violation of Bell's inequalities [7] are profoundly related 
to quantum entanglement [8]. Other applications could 
be linked directly to a resource named quantum discord 
[9-11], which is identical to entanglement for pure states 
but differs for statistical mixtures [12-14]. While for en- 
tanglement the term quantum correlation is suitable, not 
least in view of its connection to nonlocality, quantum 
discord indicates the presence of non-commuting local 
observables in the decomposition of the state which does 
not necessarily imply strong correlations [15-20]. How- 
ever, regardless of its interpretation in terms of corre- 
lations, quantum discord has proven to be an important 
resource for certain tasks in quantum communication and 
computation. It is considered especially promising in 
the context of operations involving highly mixed states, 
which emerge naturally due to the inevitable influence of 
noise [21]. 

Several methods have been developed which allow for 
the detection of quantum discord with relatively small ef- 
fort if all subsystems are under sufficient degree of control 
[15, 22-24]. Recently, we have shown that the quantum 
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discord of a bipartite system can be witnessed by access- 
ing only one of the two subsystems [25]. The method 
extends a general theoretical scheme for the detection 
of initial correlations in the dynamics of open quantum 
systems developed in Ref. [26], which has been recently 
realized experimentally [27, 28]. Typically, an open quan- 
tum system represents a well-controlled quantum system 
which is coupled to a complex, largely inaccessible en- 
vironment and therefore constitutes a natural setting in 
which we could benefit from the method described in this 
paper. 

Our strategy for the construction of a local witness 
for the quantum discord in a bipartite system is based 
on a local dephasing operation, describing measurements 
carried out on one of the subsystems, which leaves the 
marginal states invariant while erasing all quantum dis- 
cord between the two subsystems. When the subse- 
quent time evolution of the composite, bipartite system 
is changed by this dephasing operation, one can conclude 
that the original state has a non-vanishing quantum dis- 
cord. A suitable local witness for quantum discord is 
thus given by any appropriate measure for the distance 
between the time-evolved reduced subsystem states ob- 
tained from the total system states corresponding to the 
evolution with and without local dephasing operation 
[25]. 

In the present paper we develop the details of this 
method and provide a study of its applications to ther- 
mal equilibrium states of generic complex quantum sys- 
tems. In order to assess the performance of our witness 
for quantum discord we compare the actual dynamics un- 
der randomly generated Hamiltonians with the mean val- 
ues and fluctuations obtained from the average over the 
unitary group equipped with the Haar measure, employ- 
ing results of Ref. [29] . We conclude with the formulation 
of a general ergodicity-type hypothesis which relates the 
average of the local witness over the unitary group to the 
time average of the witness obtained for a generic system 
dynamics. 
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II. QUANTUM DISCORD AND LOCAL 
DEPHASING OPERATION 

Throughout this paper we deal with a bipartite Hilbert 
space H — Ha ® Hb, composed of local Hilbert spaces 
Ha and Hb with dimensions cLa and ds, respectively. 
A state p of the composite system has zero discord with 
respect to subsystem A if and only if it can be written as 
[14] 

p = J2p^^®p 1 b, (1) 

i 

with a basis of Ha, a probability distribution {pi}, 
and a set of arbitrary quantum states {pg}- States of 
zero discord are considered as classical. In the following 
we use this asymmetric definition, expressing classical- 
ity with respect to subsystem A. The reduced density 
operator pa — i s obtained from p via the par- 

tial trace over subsystem B. We introduce the quantum 
operation 

X>)(i|X A |i)<i| (2) 

i 

which represents a completely positive and trace preserv- 
ing linear map acting on operators Xa of subsystem A. 
The definition in Eq. (1) is then equivalent to the follow- 
ing statement: A state p has zero discord if and only if 
the operation 

($®I B )p = J2 n iP u i ( 3 ) 

i 

leaves the state invariant, where we have introduced the 
local projectors IT = \i)(i\®lB onto the eigenbasis of pA, 
and Ib denotes the identity operation on subsystem B. 
Equation (3) defines the local dephasing operation in the 
eigenbasis of pa- This operation constitutes the central 
element for the local detection scheme and has a series 
of important properties [25]: 

(i) The operation (3) can be interpreted as a nonse- 
lective measurement in the eigenbasis of pa, which 
is fully accessible from p by measurements in the 
local subsystem A. 

(ii) None of the two reduced density operators pa and 
Pb is affected by application of the local dephasing. 

(iii) The state produced by the local dephasing opera- 
tion is always classical. 

Property (i) is easily confirmed: Assume that the re- 
duced state pa has been obtained by state tomography. 
After diagonalization this yields the local eigenbasis {\i}}- 
The nonselective measurement in this basis is described 
by the operation (2), which by extension to the total 
Hilbert space results in the local dephasing operation (3) 
associated with the state p. Furthermore, this operation 
describes complete decoherence in the basis {\i}}- The 



diagonal elements of any operator represented in this ba- 
sis are left unchanged while all off-diagonal terms are set 
to zero. 

To prove property (ii) we use a completely general de- 
composition of p = J2 a R1 ® R% mto arbitrary fixed 
operators on Ha and Hb- The state after applica- 
tion of the local dephasing operation will be denoted by 
p' = (3>®Ib)p- Its corresponding reduced density opera- 
tor p' B of subsystem B will be unchanged, since only the 
identity operation is applied to this part of the Hilbert 
space: 

OL 

= X>{*(J*a)}J8 

a 

= Y,^{K}Rb = ^aP = Pb. (4) 

a 

The reduced state of subsystem A is not altered since the 
measurement is performed in its own eigenbasis: 

P ' A = Tr bp =^bJ2J2 ® 

a i 
a i 

= ^2Pi\i){i\= Pa, (5) 

i 

where p t = £ Q Tr{i?g} (i\RX\i) = (i\p A \i). 

Finally, property (iii) is obvious since p' can be readily 
cast into the form of Eq. (1) with pip B — 
The combination of all three properties leads to an ad- 
ditional interpretation: Performing a nonselective mea- 
surement in the local eigenbasis, i. e., applying the cor- 
responding local dephasing operation to a state p erases 
the quantum discord in p while leaving its marginals un- 
changed. 

We end this section with remarks on two special situ- 
ations. First, in the case of degeneracies in the spectrum 
of pa, the local eigenbasis is not uniquely defined. In this 
case every eigenbasis can be used for the local detection 
of quantum discord. If experimentally feasible, different 
bases can be tested to find the basis for which the wit- 
ness, to be introduced in the next section, is maximal. 
Second, if the local state tomography yields a pure state 
Pa = \ip)(ip\, no further action is required. It is already 
safe to conclude that no total correlations exist between 
the two subsystems and the total state is a product state, 
p = \(p)(ip\®pB- Specifically, this situation is encountered 
in experiments if one of the subsystems is prepared in a 
pure state. Total correlations in terms of the distance to 
the corresponding product state can be witnessed on the 
basis of an arbitrary local operation using the contraction 
property of the trace distance [26]. 
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III. LOCAL WITNESS FOR QUANTUM 
DISCORD 

Since the states p and p' differ only in their quantum 
discord, a possible measure for the amount of discord is 
given by the distance [30] 

V(p) = \\p-p>\\ 2 . (6) 

Any operator distance establishes an adequate measure. 
For later applications we choose the squared Hilbcrt- 
Schmidt norm \\A\\ 2 = TrA^A. With this choice, Eq. (6) 
can be written as a difference of purities [25] . More gener- 
ally, for any map of the form $(Xa) — J^. -KiXj^Ki with a 
complete set of mutually orthogonal projection operators 
TTi we have: 

\\p - ($ ® l B ) HI 2 = Tr {p 2 } - Tr {[($ ® I B ) p] 2 } 

= V(p)-V((<f(g>I B )p), (7) 

with the purity V(p) — Ti{p 2 }. To prove this relation 
we write the left-hand side of this equation as 

Hp - ($ ® i B ) pf 

= V{p) - 2 Tr{p($ ® l B )p} + l B )p) ■ (8) 

Making use of the Kraus representation of we obtain 

V ((* ® I S )p) = Tr J Y SijnKknRfoi ® ^< > 
= ^ Tr {j^TT,^} Tr 

= Tr J Y ® > 

= Tr{p($®I B )p}, (9) 

which proves Eq. (7). Moreover, this adds a nice oper- 
ational interpretation to the measure T>(p) in terms of 
the purity-decreasing effect of the local dephasing opera- 
tion. Note that if the state p is pure, the expression V{p) 
yields the generalized concurrence [25] , a well-known en- 
tanglement measure [31, 32], illustrating the equivalence 
of discord and entanglement in the case of pure states. 

The goal of our local detection method is to derive 
information about the total state only by means of op- 
erations on one of the subsystems, say subsystem A. To 
this end, we must take into account the dynamics of A. 
At any given moment in time, the reduced density ma- 
trix contains the full information about the correspond- 
ing subsystem. However, the same reduced density ma- 
trix can be derived from very different states of the total 
system. Consider the standard example of a maximally 
mixed state pa = 1a/ d A- This density matrix may corre- 
spond to a maximally mixed total system state p = 1/d of 
H, where d = d,Ad B . It may also stem from a maximally 



entangled pure state of H. The reduced density opera- 
tor at a single moment in time does not reveal which of 
the two cases apply. However, the nature of the total 
state can become amenable to operations carried out on 
subsystem A if we make use of the dynamics of pa- 

For simplicity, we assume that the composition of the 
systems A and B forms a closed system. We will see 
below that this assumption can be dropped. The dy- 
namics of a closed system is described by a unitary time 
evolution operator U t , propagating states from time 
to time t. Tracing over subsystem B yields the reduced 
density matrix at time t, p^(i) = Ti B {U t pU}}. The 
local detection method is based on the following idea: 
First, the accessible part of the unknown initial state 
p is measured, yielding the state pa and its eigenbasis. 
After the reference state p' is produced by local nonse- 
lective measurement of the total state in this basis, we 
compare the dynamics of the two reduced states pAif) 
and p' A {t) — Trs{?7tp'?7j }. Again, the difference of these 
states can be quantified by the distance 

dist(i) = \\ PA {t) - P ' A {t)\\ 2 = \\Tr B {U t (p P ')U}}\\ 2 . 

(10) 

First, note that dist(0) = due to property (ii) of the 
local dephasing map. On the other hand, if we find an in- 
stant of time t > for which dist(t) > 0, we can conclude 
that p and p' must be different states. This in turn indi- 
cates that T>(p) > which enables us to locally witness 
bipartite quantum discord [25]. 

We note that this method can even be extended to 
general linear time-evolutions given by a family of quan- 
tum dynamical maps A t , such that p(t) — A t (p) and 
p'(t) = At(p'), which yields 

dist(t) = ||Tr B {A i (p-p')}|| 2 - (11) 

Thereby the scheme can be used to detect correlations 
also in bipartite systems under additional dissipation 
caused by the coupling to an external environment. For 
the rest of this paper, we will restrict to the case of uni- 
tary evolution. 



IV. PERFORMANCE OF THE WITNESS AND 
EXAMPLES 

The above method may fail to detect correlations de- 
pending on the time evolution U t . Consider for instance 
the trivial case of two uncoupled subsystems. The time 
evolution factorizes, U = Ua <S> U b , where we omit the 
time argument. In this case, no signature of the to- 
tal state will be visible in the reduced system dynam- 
ics, which can be seen easily by decomposing p — p' — 
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Tr B {U(p - p')U^} = Tr s |E U A°% U A ® U BD%Uy> 
= Y,UaD%U]^{D%} 

a. 

= C/ A Tr B {p - p'} U\ = 0. (12) 
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The question is thus, what is the performance of the 
method for generic systems? In order to answer this ques- 
tion we make use of a recently developed approach based 
on unitary average values [25, 29]. In order to obtain an 
estimate for the quantity dist(i), we replace Ut with a 
random unitary matrix U and determine the average in- 
tegrating over the uniform Haar measure dp. According 
to ensemble theory, the average value is expected to re- 
flect the behavior of generic complex quantum systems. 
We denote unitary average values by angular brackets, 



(F(U)) = / dp(U)F{U). 



(13) 



The Hilbert-Schmidt distance for an arbitrary pair of 
states p and p' yields the average value [25] 



\Tr B {U(p-p')tf}\ 
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FIG. 1. (Color online) The plot shows the dependence of the 
unitary average value fi and the variance s 2 on the parameter 
z for the Werner states. The relative error is constant at 
s/m~ 0.58. 



A. Werner states 

As a first simple illustration of this method, we con- 
sider the Werner states. 



(17) 



with |*) = (| 00) + |11))/V2 and -1/3 < z < 1. The 
purity yields V(pw) = j(l + 3z 2 ). The reduced system 

state is the maximally mixed state pa = 5(|0)(0| + |1)(1|), 
and to produce the reference state by local dephasing we 
can project onto the operators |0)(0| and |1)(1|, 



(M) f/ = ($®I B )p = 53(|i)<i|®lB)p(|i)<i|«8lB). (18) 



and the variance [29] 



Var(||Tr B {U(p p')tf}\\ 2 ) = Cl (Tr{(p - p'f}f 

+ C2 Tr{(p-p') 4 }, (15) 

with the coefficients c\ and ci given by 



(-2 



444 + dU%)(d 2 



2(15 -™ A » B 

(36 - lM\d B + d\d%)(d 2 A d% 
-10d A d B (d B -l)(d A -l) 



I) 2 



44(44-7)2-36 



(16) 



Inserting p' = ($ ® into Eq. (14), we find that 

the average increase of the local distance is directly pro- 
portional to the squared Hilbert-Schmidt distance of the 
original state p to its locally dephased reference state p' , 
which we had previously defined as T>(p), a measure for 
quantum discord. This result also holds for a more gen- 
eral average, which is performed only over the eigenvec- 
tors of the Hamiltonian while the time dependence and 
the eigenvalue distribution are retained, see Refs. [25, 29]. 
Hence, we conclude that for generic open quantum sys- 
tems the quantum discord in the initial state will be suc- 
cessfully detected by the present method [25]. 



Thus, we obtain the reference state 
1 + z 



($®l B )p w 



(|00)(00| + |11)(11|) 



1 - z 



(|01)(01| + |10)(10|), (19) 



with the purity V((<$>®I B )pw) = \{^ + z 2 ). The generic 
increase of the distance in the reduced system is given by 
Eq. (14), which leads to 



p = (||Tr B {U(p w - (3>®I B )p w )U*}\\ 2 ) = l -z- 



(20) 



where we have made use of Eq. (7). 

The Werner state has nonzero discord for any value 
of z 7^ 0. On the other hand, the state is separable 
for z < 1/3, which illustrates that quantumness beyond 
the concept of entanglement can be identified with this 
method. The variance is given by Eq. (15), which for this 
state yields 



Var 



(||Tr B {C/(p-(*®lB)p)f/ t }| 



19 
1400 



(21) 



The relative error is constant for all values of z and 
amounts to s/p — \/l9/56 ~ 0.58. The relatively large 
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FIG. 2. (Color online) Comparison of the unitary average 
with the actual time evolution for the Gibbs states of six 
randomly picked two-qubit Hamiltonians {cLa — ds = 2) at 
fixed temperature {3=1. The pictures show the value of the 
Hilbert-Schmidt distance after applying the local detection 
method to the Gibbs state. 



value of the variance is explained by the low dimensions 
of system and environment. A plot showing the depen- 
dence of expectation value and variance on the parameter 
z is given in Fie. 1. 



Random Gibbs states of 2 x do systems 



In this section we demonstrate the local detection 
scheme for Gibbs states of randomly generated d- 
dimensional Hamiltonians. Once such a random H has 
been generated [33], the Gibbs state can easily be ob- 



tained as pa 



/Z, with the partition function 



Z = Tre _/3H , /3 = 1/kT, temperature T, and the Boltz- 
mann constant k. We consider the total Hilbert space 
to be 2rfs-dimensional, i. e., the system Hilbert space 
% A is two-dimensional. Employing the product basis 
{|0) , |1)} ® {|Xi)},=ij where {|xi)} denotes an arbitrary 
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FIG. 3. (Color online) Comparison of the unitary average 
with the actual time evolution for the Gibbs states of four 
randomly picked Hamiltonians for a qubit coupled to environ- 
ments with different dimensions at fixed temperature {3 = 1. 



fixed basis oIT-Lb, the Gibbs state pa can be written as 

PG = J>™ |o)(o| ® \ Xi ){ Xj \ loX 1 ! ® \xi){xj\ 

(22) 



'•J 



>■:) 



Hence, the reduced density operator of subsystem A can 
be represented by the matrix 



p A = TtbPg 



(23) 



On the basis of the eigenvectors {|0) , |1)} of this (2 x 2)- 
matrix, the local dephasing map is expressed as 



(* ® i b )p = n^n^ + n T pn T , 



(24) 



with lb = \i)(i\®1b- Application of this map to the origi- 
nal Gibbs state pq creates the reference state (<J?®Ib)pg- 
Next, we examine the dynamics of the distance of the two 
reduced system states by creating the corresponding time 
evolution operator Ut = exp{—iHt} from the same ran- 
domly generated Hamiltonian H. The distance is given 
as a function of t by: 



dist(i) - \\Tr B {U t (p G - (^®Ib)pg)UJ}\ 



(25) 



On the other hand we can obtain the unitary expectation 
value and its variance for the same quantity by Eqs. (14) 
and (15), which in this case yield 



(i = (\\Tr B {U( PG - ($g> I b )pg)U^}\\ 2 ) 
3d B 

T 



44 



\p G - ($® I B )pc 



(26) 
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FIG. 4. (Color online) Dependence of the discord of a ran- 
domly picked fixed thermal state on the temperature. For 
higher temperatures (lower /?), discord diminishes. The pic- 
ture on the lower right shows the average value and one stan- 
dard deviation as a function of the inverse temperature /3. 



and 



s 2 = Var (||Tr B {U(p G - (3> ® I B )p G )f/ t }| 

3(15-16d| + 164) „ ,^ „, l|4 
^ 2(l-44) 2 (44-9) IIPG ' ( ^ 1Ii3)PGl1 

(27) 

Figure 2 shows a series of time evolutions including 
the corresponding unitary average value /i and the first 
standard deviation s for six randomly generated 2x2 
Hamiltonians at fixed temperature j3 = 1. The depen- 
dence on the environmental dimension is plotted in Fig. 3, 
while Fig. 4 displays the role of the inverse temperature 
/?. From these simulations we can make a number of ob- 
servations. First, we see that generic Gibbs states con- 
tain quantum discord since the function dist(i) assumes 
nonzero values for all realizations, confirming measure- 
theoretic studies on the abundance of quantum discord 
[20, 34]. Second, for most of the examples the time evo- 
lution fits nicely into the margin given by the unitary 
average within one standard deviation, indicated by the 



highlighted areas. It is of course no surprise to find some 
deviating realizations as in the top right picture of Fig. 2. 
Third, as becomes obvious by comparison of Figs. 2 and 
3, the unitary average value depends stronger on the di- 
mensions of system and environment than on the actual 
Hamiltonian. The values in Fig. 2 differ only very lit- 
tle between the considered random examples, while in 
Fig. 3 we see that the average value [i and the stan- 
dard deviation s decrease significantly with increasing 
environmental dimension. This is mainly caused by the 
dimension-dependent factors in Eqs. (26) and (27). 

Figure 4 shows how the witness dist(t) changes for dif- 
ferent temperatures. We see however that the overall 
functional shape remains similar which is due to the fact 
that the Hamiltonian is the same in all plots. The bot- 
tom right picture shows the asymptotic convergence of 
mean value and variance for decreasing temperature. In 
the high-temperature limit (/? — > 0) the unitary average 
value, and with it the generic effect of the initial corre- 
lations on the reduced system vanishes as expected since 
the state becomes closer to a complete mixture, which 
is a state of zero discord. Note that correlations in the 
low-temperature limit of the Gibbs state can be used to 
reveal the structure of the ground state [35], which in 
turn can be associated with a quantum phase transition 
[36-39]. 

To conclude this section, we recall that a state of 
nonzero discord cannot be a factorized product state [40] . 
On the other hand, factorizing initial conditions are com- 
monly assumed in the derivation of master equations for 
the dynamical description of open systems in terms of 
completely positive maps, see, e.g., Refs. [21, 26, 41] and 
references therein. Hence, the present method can also 
be used to detect deviations from this assumption [25]. 
Obviously, if the witness is nonzero, a dynamical map 
which is independent of the correlations does not exist. 
A study of the role of the total initial correlations in 
thermal equilibrium states is presented in Ref. [35]. 



C. An ergodicity-like relation 

The foregoing study shows that unitary averages pro- 
vide important and useful information about the time 
evolution, which may be experimentally observable. It 
was pointed out in Ref. [29] that the dimension ds ap- 
pearing in expressions for the averages must be chosen 
carefully. Formally, it is always possible to artificially in- 
crease the dimension of the Hilbert space by including an 
additional Hilbert space which is not coupled to the orig- 
inal system. Correspondingly, the dimension appearing 
in the expectation value must be regarded as an effec- 
tive dimension, indicating the dimension of the subspace 
of the Hilbert space which actually affects the local dy- 
namics. Thus we are led to an ergodicity-like hypothesis 
about the equivalence of the unitary average value ob- 
tained with a suitable, effective dimension d c g and the 
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time average according to the given, actual Hamiltonian: 



\TT B {U(p-p')U^\ 



cir 



lim — 

T-yoo T 



dt 



Tr, 



{u t (p-p')U}}f . (28) 



For complex generic systems, the effective dimension co- 
incides with the dimension of the Hilbert space. The 
effective dimensions depend not only on the system pa- 
rameters but also on the observable in question. For 
example, in a partly chaotic system with regular areas, 
some initial states may explore large parts of the state 
space in the course of their time evolution while for dif- 
ferent initial conditions only a very limited fraction may 
be visited. The estimation of the dimension of quantum 
systems is a topic of growing interest [42]. 



systems and quantum communication protocols. The 
procedure was illustrated by application to thermal equi- 
librium states of random Hamiltonians. In order to esti- 
mate the performance of the method for generic systems 
we compared the time evolution with the value obtained 
by averaging over all unitary evolutions employing the 
Haar measure. The mean values as well as the fluctu- 
ations predicted by the Haar measure were found to be 
in good agreement with the actual time evolution. This 
fact led to the proposition of an ergodicity-like hypoth- 
esis, linking unitary average and time average via the 
introduction of an effective dimension of the underlying 
Hilbert space. Further studies are required, on the one 
hand to confirm this hypothesis with additional examples 
of generic systems and, on the other hand, to obtain the 
effective dimensions of non-generic systems which typi- 
cally exploit only an effective subspace whose dimension 
is much lower than that of the total Hilbert space. 



V. CONCLUSION 

The method discussed in this paper allows for the de- 
tection of quantum discord in bipartite systems when ac- 
cess to only one of the subsystems is possible. This situ- 
ation emerges naturally in the context of open quantum 
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